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Abstract 

We report on an investigation of the LANL method for determining the 
0(a) improvement coefficient ca nonperturbatively. We find we are able to 
extract reliable estimates for the coefficient using this method. However, our 
study of systematic errors shows that for very accurate determinations of ca, 
the smearing function must be tuned and the volume fixed to keep the 0(a) 
ambiguity in ca fixed as /3 varies. Consistency was found with previous results 
from the LANL group and (within fairly large errors) 1-loop perturbation 
theory; ca does not change significantly over the range (3 = 5.93—6.2. The 
big difference between our results and those of the ALPHA collaboration, 
around [3 = 6.0, show that the 0{a) differences in ca between the different 
methods can be large. We find that the lattice spacing dependence of f^^ and 
the renormalised quark mass is much smaller using our values of the coefficient 
compared to those of the ALPHA collaboration. 

I. INTRODUCTION 

The use of Symanzik improvement |l| of lattice actions and matrix elements is widespread 
and very effective. However, with each improvement term added the corresponding coeffi- 
cient must be determined to enable discretisation effects to be reduced to the desired level. 
Considering the light hadron spectrum and matrix elements, the relevant 0(a) improvement 
coefficients are, for the most part, only known to 1-loop in perturbation theory, leaving 
residual 0{a^a) discretisation terms. A nonperturbative determination of these coefficients 
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is desirable to completely remove 0{a) effects. Such a determination is possible through the 
imposition of the axial Ward identities (AWI) on the lattice. 

Central to a programme of determining 0{a) improvement coefficients nonperturbatively 
is the calculation of ca, the improvement coefficient of the axial-vector current; the improved 
current appears in the expression for the generic axial Ward identity and ca must be de- 
termined before other operator improvement coefficients can be calculated 0. So far, two 
groups, the LANL group I^J^Q and the ALPHA collaboration 0,^, have calculated ca along 
with several other improvement coefficients. 

Their results for at /3 = 6.0 and 6.2 are summarised in figure |I| and compared with 
1-loop perturbation theory (using a = ap{l/a) [0]). While the results are compatible on 
the finer lattice, there is a big difference in the values for (3 = 6.0. This difference can be 
explained by the 0{a) ambiguity which exists in nonperturbative determinations of 0{a) 
improvement coefficients. As long as the improvement conditions for each determination 
are applied consistently as f3 changes, differences in the value of ca are not important in 
principle; the difference will disappear in the continuum limit. 

In practice, different values of ca can have a large effect away from the continuum limit. 
This is because the matrix elements appearing at 0{a) for and the renormalised quark 
mass are numerically large compared to the leading order term. 

nP={P\A,\0) + CA{P\ad,P\0) (1) 
_ {P\d,A,\0) , {P\adlP\0) 
^PCAc - 2{P\P\0) + 2{P\P\0) ' 

and a {P\d4P\0) / {P\A4\0) ~ aM^/niq. An 0(a) ambiguity in ca therefore appears at O(a^) 
but multiplied by a large matrix element and it is undesirable to have large O(a^) scaling 
violations, even if 0(a) errors have been removed. 

Our aim is to investigate how well ca is determined using the LANL method. The 
latter only requires a conventional analysis, which is available from simulations performed 
for spectroscopic calculations, compared to using the Schrodinger functional techniques of 
the ALPHA collaboration. With significantly higher statistics than those of reference 
we are able to improve on the LANL analysis by performing correlated fits, investigating 
the choice of lattice derivatives employed more widely and determining the stability of the 
LANL results to changes in the fitting range. In addition, we study the scaling behaviour 
of and the renormalised quark mass with respect to the choice made for ca- 

The paper is organised as follows: in section ^ we outline how to extract ca from the 
PCAC relation and in particular the method employed by the LANL group. Results are 
presented in section |ITT|, which includes a comparison of our results with those of the LANL 
group and the ALPHA collaboration. The scaling of and the renormalised quark mass 
is dealt with in section followed by the conclusions in section 0. Technical details not 



directly related to the method of calculating ca - the simulation details, extracting meson 
masses and decay constants, the renormalisation factors used to obtain f^^ and m*^"^ and the 
chiral extrapolations are all given in the Appendix. 
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II. Ca from the pcac relation 



The PCAC relation in euclidean space can be written as 

{Jd^A^{x)) = 2mpcAc {JP{x)) (3) 

and should hold on the lattice for all x not coincident with J up to discretisation terms. The 
axial- vector current, A^, is given by ip'j^'j^ip, the pseudoscalar operator P = 'ip'-y^'ip and J 
is any operator with the pseudoscalar quantum numbers. mpcAc is the bare current quark 
mass. For simplicity we sum over position space, restricting ourselves to zero momentum, 
and define 

<JdMtl>_ 
"^'^^^ <JP{t)> • 

Thus, equation || becomes 

rj{t) = m{t) = 2mpcAC + 0{a). (5) 

This relation holds for all states (ground and radial excitations) of the pseudoscalar meson. 
In the limit of large times, when only the ground state contributes to rj, then m{t) is a 
constant given by 2mpcAc + O^.^Xo), where Og.s.(a) are the discretisation errors associated 
with the ground state. At earlier times, when excited states make a significant contribu- 
tion (with different discretisation errors), m{t) becomes time dependent, as can be seen in 
figure |. 

The size of the discretisation terms (and the time dependence of m{t)) are reduced to 
O(a^) when we improve the axial- vector current [| 

A^^ Ai = A^ + acAd^P + 0{a'^). (6) 

Then 

rj{t) + acASj{t) = niimpit) = 2mpcAc + 0{a^), (7) 

with 

< JdlPjt) > 
''^'^ = < JP{t) > ■ 

Clearly, changing J or the time t changes the size of the contribution of each state to rj, 
and the size of 0(a) terms. The improvement term must still cancel these terms, however, 
giving rise to a quark mass which differs from rriimp only in O(a^). 

rj(t') + acASjit') = niimpit') = 2mpcAc + O'(a^), (9) 



^The quark action must also be improved to O(a^) using the Sheikholeslami-Wohlert term with 
the value of cgw determined nonperturbatively. 
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or 



rj>{t) + acASj'{t) = m-„p(t) = 2mpcAC + 0"{a ). (10) 
By forcing niimpit) and rriimpit') to be equal we can solve for ca- 

1 r(t') - r(t) 



as{t') - s{t) 



CA. (11) 



A suitable choice for t and t' is to set t = tgs in the region where the ground state dominates 
rj and sj and t' = tex in the region where there is significant contribution from excited 
states. 

In order to illuminate the 0{a) ambiguity in ca, we identify 

Arj = r,j{t') - r,j{t) = Ao[0(a)] + Arj[0(a2)] + . . . (12) 

Asj = sj{t') - sj{t) ^ -Arj[0{a)]/acA (13) 
l rj{t')-rj{t) ^ Ao[0(a^)] 

a5,(t')-3j(t) ^ ^ ^Ao[0(a)] ^'^^ 

where Arj[0(a)] denotes the change in rj due to 0(a) effects etc. Thus, the 0(a) ambiguity 
in Ca depends on the difference of the O(a^) (and 0(a)) terms between r at tex and tgs rather 
than the absolute values. Obviously if Arj[0(a^)] ~ Arj[0(a)], the error in ca will be as 
large as ca itself (and the above expansion will not be valid). 

The LANL method is equivalent to using equation |TI|. It involves performing a fit to 
r.]{t) and s.j{t) such that rj{t) + acASjit) is equal to a constant (2m), where ca and 2m are 
parameters in the fit. The fitting range is chosen to be from tgx to tgg. The advantage of 
performing a fit over calculating the ratio in equation ^ is that one can test the ansatz that 
the value of ca reduces the time dependence of rj{t) (and hence the discretisation errors) 
with the x^. 

We note that the ALPHA collaboration employs a slightly different method to calculate 
Ca- Within the Schrodinger functional approach it is possible to simply work in the region 
when the ground state dominates and rj and sj have plateaued. The boundary fields are 
varied and this changes the discretisation errors in the ground state. A ratio similar to 



equation [TT| is built up from rj and sj from different boundary fields. A feature of using 
the Schrodinger functional technique is that the analysis can be performed at directly zero 
quark mass. 

The unknown 0(a) ambiguity in ca is due to 0{a?) errors in the axial-vector current, 
the pseudoscalar current and the light quark action. The improvement scheme which we 
have chosen defines ca in the limit of zero quark mass, however, ca is calculated at finite 
TUq and then extrapolated. There is an additional 0(a) ambiguity due to the discretisation 
chosen for the temporal derivatives in equations ^ and H, which is proportional to niq and 
vanishes in the chiral limit. As part of our analysis we investigated two different choices for 
the lattice derivatives in the determination of ca- "standard" symmetric lattice derivatives 

d, ^ A(+-) = i (5,v^+A - (15) 

A^ ^ 5x,x-\-ji '^^x,x ^x,x—ji 
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which contain O(a^) errors, and "improved" 0{a'^) derivatives 



A|+-) = A(+-) - iA+A(+-)A- (16) 



A(2) = A(2) - 1 
M 12 



A+A- 
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where A"*" = — 1 and A^ = 1 — Two points should be taken into account when 

choosing the form for the derivatives. As one improves the derivatives then the minimum 
value that t^x can take becomes larger since there must be no overlap with the origin. In 
addition, noise rapidly dominates the determination of ca as tex increases and there is only 
a small window of timeslices from which to extract ca- The LANL group considered another 
form of O(a^) lattice derivative which has smaller O(a^) terms compared to eqns ^ 



but uses fewer timeslices than the O(a^) derivatives. While the different definitions give 
consistent results for ca in the chiral limit, we found using improved derivatives helped in 
extracting ca- This point is discussed in the next section. 



In addition to the LANL method and equation [IT] we considered extracting ca by chang- 
ing J from ipj^tp to V'7475V'- The change in J does not give rise to significantly different 
discretisation effects and the determination of ca was not improved. In the following we set 
J = P and drop the subscript on r and s. We also used the LANL method with finite mo- 
mentum correlators in r and s (including the additional spatial derivative terms). However, 
at finite momentum the 0(a) errors in ca are increased and no additional constraint on the 
coefficient is obtained compared to the zero momentum results. 



III. RESULTS FOR Ca 

We have tested the LANL method using the UKQCD quenched data set. The simulation 
details are given in table | and discussed in the Appendix. The best analysis was possible 
at p = 5.93. 



A. Results at /? = 5.93 

To implement the LANL method we first fix tgg. Figure ^ shows the fractional contri- 
bution of the ground state to the correlators which appear in r and s with ki = 0.1327 
for the two types of correlators, LL and FL, that are available for this data set (see Ap- 
pendix P). The ground state dominates by approximately timeslice 12 for all correlators. 
We use tgs = 12 and as a check also tgs = 14. tex is allowed to vary in the region tgx < tgs 
while being careful to avoid any overlap with the origin. 

We present the details of the calculation of ca for the heaviest k = 0.1327 using standard, 
O(a^), derivatives in table |I|. The differences, s(tex) — s{tgs) and r{tex) — r{tgs), are well 
determined for t^x close to the origin, but rapidly decrease and become dominated by noise 



as tex increases to around timeslice 8. The ratio of equation |TT| agrees with the results 
obtained using a fit apart from tex = 3. However, the LANL fit does not give a reasonable 
Q (defined as Q > 0.01) until timeslice 5 for FL correlators and 7 for LL correlators. 
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The LANL method is expected to achieve good fits in a region where the O(a^) errors 
in rriimpit) (equation 0) are not rapidly changing. In this region, the time dependence of 
acAs{t) can compensate for any corresponding variation in r(t). This is hkely to occur when 
there are only a few states contributing to r and s 0. Ideally, in this region, the estimates 
of ca are stable with a variation of tex- Any dependence on tex, or any difference between 
the results from the LL or FL correlators means that O(a^) contributions to r and s are 
appearing in ca- 

The results from the LANL fits, with reasonable Qs, are shown in figure |[ ca is stable 
with tex (possibly excluding t^x = 5 for FL correlators), although noise rapidly dominates 
as the fitting range becomes smaller. There is agreement between the values obtained using 
FL and LL correlators, and tgs = 12 and 14. We can compare the range over which a 
good fit is found with the fractional contribution of the sum of the ground and first excited 
state to the correlators which make up r and s, shown in figure |^. Roughly, the earliest 
tex with Q > 0.01 corresponds to the timeslice when all but the first excited state and the 
ground state dominate the correlators; since the FL correlators have a lower contribution 
from excited states, the ground state plus first excited state dominate at an earlier timeslice 
compared to the LL correlators and a smaller t^x can be used. 

The results change quite significantly when we switch to improved derivatives. Figure |^ 
shows the effects of changing the derivatives on m(t) = r{t) and s{t) for LL correlators. 
Clearly, the time dependence of m{t) is much reduced in the range t = 5 — 10 when improved 
lattice derivatives are used, indicating that most of the discretisation effects seen when using 
the standard derivatives are due to the O(a^) error in r{t) associated with this derivative 
and not the 0{a) errors which we are trying to cancel with acAs{t). A similar but much less 
dramatic effect is seen for s(t). This translates into much smaller values for r(tex) — r{tgs), 
s(tex) — s(tgs) and ca, as seen in table fT| and figure ^. 

For both LL and FL correlators, reasonable fits can be obtained with slightly smaller 
tex than with standard derivatives, suggesting Arj[0{a^)] in this region is not large. As 
figure ^ shows, there is agreement to within 2a between the LL and FL results, with the 
exception of the fit to LL correlators with 5 — 12, which disagrees significantly with the FL 
result over the same fitting range. This is presumably due to a larger 0(a) in ca for the 
LL result, since these correlators have a larger contribution from first excited (and higher) 
states. This fit is on the borderline of being considered reasonable; changing tgs to 14, the 
Q drops further. 

We demonstrate the effect of 0(a) improvement on 2ampcAc, for various values of ca 
for LL correlators in figure |^. The discretisation errors in m{t) can be reduced using either 
standard or improved lattice derivatives, however, the latter requires a smaller value of ca- 
As expected rriimpit) is constant over the time range used in our fit, with the plateau being 
one timeslice longer in the case of improved derivatives. 

Our improvement condition is defined in the chiral limit and the results for ca must be 
extrapolated to zero quark mass. Details of this procedure are given in Appendix and the 



^As (3 increases and the O(a^) errors decrease one expects a reasonable fit to be possible including 
more excited states than on coarser lattices 
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results are presented in figure ^ and in table |^ In general, ca from using standard deriva- 
tives has a bigger statistical error than that obtained using improved derivatives since t^x 
is larger and the chiral extrapolation is usually more severe. The latter point is illustrated 
in figure ^ From the discussion above, the large value of ca from standard derivatives is 
due to a large contribution to rriimpit) from the O(a^) errors associated with these deriva- 
tives. These errors are rUg dependent and should disappear in the chiral limit. The figure 
shows that ca drops significantly with quark mass, agreeing with the result from the im- 
proved derivatives in the chiral limit. In contrast ca from improved derivatives is much less 
dependent on the quark mass. 

Comparing all results in figure ^ we find ca{k,c) from different derivatives, smearing and 
from different t^x and tgs are consistent in the chiral limit (with the exception of tgx = 5 for 
FL correlators) . We also implemented the LANL lattice derivatives and found consistent 
results in the chiral limit, ca = —0.032(14) from LL correlators, fitting 6—14, using improved 
derivatives is taken as the final value for ca at this (3. The error reflects the spread in values 
for Ca and indicates the uncertainty from some of the associated 0{a) effects. 

We now consider applying the LANL method at different /5s. In principle, one should 
ensure, as accurately as possible, that the same improvement conditions are applied to 
determine ca as (3 is changed. This ensures the systematic errors are correlated between 
different determinations and ca smoothly extrapolates to zero in the continuum limit. For 
example, we need to keep the proportion of the excited states to ground state contributing to 
r(tex) fixed. This requires a tuning of the smearing function 0, which was not possible in this 
study (and was not attempted in reference 0). Thus, we have chosen a fairly conservative 
error for ca to take into account the difficulty in applying the same improvement conditions 
for the other simulations. In principle any value of ca in figure |] is a valid estimate of the 
coefficient for a particular simulation. A more aggressive choice for ca would be, for example 
—0.050(3) from FL correlators with tex = 4. In the following, to keep the systematic errors 
as correlated as possible we extract final numbers for ca from LL correlators and improved 
lattice derivatives (as used for our choice of ca above). 

In addition, the physical volume of the simulation should also be kept fixed when deter- 
mining Ca- We comment on this in the next section. 



B. Results at /3 = 6.0 and /3 = 6.2 

Considering the analysis at /3 = 6.0 first, we present the results for ca in figure ^ and 
table 0. In addition to the simulations on the 16^ volume with FL and LL correlators, 
Ca was also calculated on a small ensemble of large volume (32^) configurations with SL 
smearing (unfortunately no LL correlators were available). As discussed in Appendix 0, 
the fuzzed smearing was optimised for the ground state and hence the first excited state 
amplitude for the FL correlators is very small; an estimate for ca could only be extracted 



^It is advantageous to work in the regime where the first excited state is the only significant radial 
excitation. A smearing with a good overlap with this state (likewise for another smearing with the 
ground state) would enable this proportion to be fixed accurately as /3 changes. 
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for one value of tex- Nevertheless, the FL and LL results are consistent and there is no 
significant variation with tgs- 

A discrepancy was found, however, when comparing to the SL results on the larger 
volume. This can be seen in figure ^ for ca_ at finite rrig around the strange quark mass. 
There is a 3.5cr disagreement between the FL value for tex = 6 and the SL value at t^x = 
4 (tg,s. = 16). The discrepancy between the results could be due to 0(a) terms arising 
from the use of different smearings or it may indicate a finite volume effect. The 16^ lattice 
corresponds to a physical volume of approximately (1.5fm)^, which is probably too small to 
accommodate excited pseudoscalar states, ca itself is an ultra-violet quantity but may be 
affected by finite volume effects because of the matrix elements being used to determine it. 

We attempted to investigate finite volume effects by comparing masses and decay con- 
stants (the matrix element < Pd^A^ >, in r(t), is related to fps) on the two volumes. 
Unfortunately, we were only able to extract these quantities for the ground state on the 
large volume. We found an .8% or 2.5(T decrease in mps changing from the small to large 
volume and no significant change in fps (this is in agreement with the results in refer- 
ence [H). A more thorough investigation of finite volume effects is needed. If the physical 
volume was the same as at /3 = 5.93, it would not matter how dependent ca is on the size 
of the lattice, since it is a higher order effect. However, the finer lattice is 16% smaller than 
that at /3 = 5.93. This must be considered when quoting an error on ca- 

The chiral extrapolations of at /3 = 6.0 proved difficult for ca from the LL and FL 
correlators. Apart from the LL result for tex — tgs = Q — 16, the errors on the extrapolated 
values are very large due to having to use a fit function quadratic in a?Mpg and/or only 
being able to fit to the lightest data points. The discrepancy of the 6 — 16 result with the 
4 — 16 SL value is approximately 2.5(T. 

As noted in the previous section, to keep the same improvement conditions ca should 
be extracted using a tex with the same relative proportion of ground state to excited states 
as that used at /3 = 5.93. One possibility, concentrating on LL correlators, is to fix tex to 
correspond to the same physical time; tex = 6 chosen at /3 = 5.93 corresponds to approx- 
imately timeslice 7 for j3 = 6.0. The statistical errors of ca{i^c) for this tex are too large 
for the estimate to be useful in our later analysis. If we choose tex = 6 then the errors do 
not reflect the unresolved 0{a) or finite volume effects mentioned above. These problems 
motivate us to discard the results for ca at this p. 

At P = 6.2 the situation is more straightforward as displayed in figure ^ and table [VI . 
There is consistency between the results from different smearings, where the fuzzed smearing 
is optimised in a similar way to that at /? = 6.0, and also as tex is varied. We also found no 
change in the results if tgg is increased and there was no difficulty with chiral extrapolation. 
Keeping the same physical tex as at /5 = 5.93 corresponds to using timeslice 9. Unfortunately, 
the estimates of ca have fallen into noise at this point. However, given the stability of the 
results with tex we are unlikely to introduce significant systematic errors if we choose a fitting 
range of 7 — 16. Thus, our final result at this (3 is ca = —0.031(5). The lattice volumes 
aX (3 = 5.93 and 6.2 are fairly close in physical size, and we assume that the error on ca is 
sufficient to compensate for the small discrepancy. 
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C. Comparison with previous results 



Figure |I] compares our results with those of the ALPHA collaboration ||^ and the LANL 
group 0,0] in the range of /3s we have simulated. Unfortunately, our errors on ca are quite 
large after chiral extrapolation. Nevertheless, we obtained consistency with the LANL re- 
sults, in particular at coarser lattice spacings. The LANL results are split up into those 
extracted using standard O(a^) derivatives and those obtained using modified {a?) deriva- 
tives, mentioned in section The latter is closer to the choice of derivatives employed 
here and we find greater consistency with our results, aX (3 = 6.2, in this case. The LANL 
results have smaller errors compared to our values even though our study has much higher 
statistics. We believe this is due in part to our more conservative error estimates to take into 
account the difficulty in applying the improvement conditions consistently as (3 changes. In 
addition, the LANL group employed much heavier light quark masses, over a wider range, 
than in our analysis, and this led to lower statistical errors after chiral extrapolation. 

The results from the LANL method are very slowly varying with (3 and do not change 
significantly from (3 ~ 6 to 6.2. This is in agreement with the 1-loop perturbative results, 
also shown in the figure for a = ap{l/a) [|^, where ca = —0.0952a QJ^l- results, 
with large errors, are consistent with the perturbative result, however the LANL values 
are somewhat higher. The uncertainty in the perturbative value, from higher order terms, 
is difficult to estimate. One could take anything between the square of the 1-loop term, 
6ca ~ 0.0006 to la^ ~ 0.04 - 0.08 in the range oi (3 = 6.2 to 5.93. The 2-loop contribution 
would have to be quite large to obtain consistency with the LANL results 0]. However, 
the calculation of this contribution would significantly reduce the uncertainty on ca] the 
perturbative result is valid in the infinite volume limit and there is no 0{a) ambiguity in 
Ca Q which is present in the nonperturbative determination and can be large, in particular 
at coarser lattice spacings. 

This can be seen when comparing our results (and those of the LANL group) with 
those of the ALPHA collaboration. At (3 = 6.2, all results are consistent, however, at 
coarser lattice spacings a large discrepancy appears as the ALPHA ca rapidly increases 0. 
This discrepancy indicates how large the 0(a) ambiguity in ca can be. In addition, the 
LANL group using new results at /? = 6.4 have found that the difference between their 
results and that of the ALPHA collaboration requires O(a^) terms as well as 0{a) [Q. We 
believe that by extracting ca in a region where only the first excited state and ground 
state contributes, looking for consistency between results from different smearings and using 
improved derivatives we have minimized the 0{a) artifacts in ca within the LANL method. 
Nonetheless, if the improvement conditions are kept fixed accurately as (3 is changed, large 
variations in the estimates of ca are not important. However, practical difficulties arise if the 
choice taken leads to significantly worse scaling violations for physical predictions, namely. 



^Although, of course, 0(a"a) terms remain in the axial-vector current. 

^At P = 5.93, CA = —0.16 from the Pade expansion of the ALPHA results j^. However, it is more 
likely to be around —0.11 @. 
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/,r and the renormalised quark mass. 

IV. SCALING OF AND 

Ca is needed, for 0(a) improved estimates of the pseudoscalar decay constant and the 
quark mass determined from the bare PCAC quark mass. Figure |T^ displays our results for 
the renormalised decay constant in units of tq as a function of the squared meson mass using 
our values for c^, at /? = 5.93 and 6.2 and also using ca determined by the ALPHA collab- 
oration at /3 = 6.0 and 6.2. The extraction of the decay constant and the renormalisation 
factors used are detailed in Appendices ^ and y. 

The figure clearly shows that using our smaller values of the improvement coefficient 
there are no significant scaling violations between jS = 5.93 and (3 = 6.2, in contrast to 
the significant violations found using the ALPHA values for ca- In the same figure we plot 
ro/^^" as a function of {a/roY for the reference mass {rQMps)"^ = 3.0. To see whether both 
data sets are consistent with the same continuum limit we fit the data simultaneously with 
a function linear in (a/ro)^ and obtain a /d.o.f. = 2.4. This is on the borderline of being 
considered a good fit. Results at more /? values are needed in order to be able to include 
higher order terms in the continuum extrapolation and accommodate the results obtained 
using the ALPHA values of ca- 

We calculated the renormalised quark mass from the improved bare PCAC quark mass. 
Details of this calculation and a consistency check, extrapolating [roMps)'^ as a function 
of fQTnpcACi given in Appendices ^ and |C[ We note that in order for the meson mass 
and the PCAC mass to vanish at the same point, chiral log terms had to be included in the 
chiral extrapolation. Similarly, we found that when these terms are included to extract 
from Mp5, consistency is obtained with extracted from rnpcAc- 

Figure [T^ presents our results for the renormalised quark mass in the MS scheme at 
the scale 2 GeV as a function of {a/roY for the reference mass [r^MpsY = 3.0. In (a) the 
calculation proceeds via the renormalisation-group invariant mass using nonperturbative 
renormalisation []lT| and the overall statistical and systematic errors are small. Again, the 



scaling violations are small using our values of c^, in contrast to the severe lattice spacing 
dependence seen using the ALPHA values. A combined linear fit to both data sets has a 
rather poor /d.o.f. = 3.3. As for the decay constant, more points are needed to check for 
a consistent continuum limit. 

In (b) we compare m^^^ determined from Kc (which is independent of ca) with that 
obtained from the PCAC mass. 1-loop perturbation theory is used for the renormalisation 
factors and hence the overall errors are much larger than in (a). The best scaling behaviour 
is seen for m^^f, followed by rrip^^^ determined using our values of c^i. A linear fit to the 
combined data gives, x^/d-o.f. = 0.1, but the errors are rather large. 

V. CONCLUSIONS 

We have undertaken a study of the difficulties and systematic errors inherent in a non- 
perturbative determination of ca using the LANL method. Some of these points may also 
apply to determinations using the ALPHA method. 
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For practical purposes it is desirable that the determination of ca is defined so that 
the improvement term for the axial- vector current removes only 0(a) errors and does not 
unnecessarily add large additional O(a^) artefacts. We conclude that this is possible with 
the LANL method but careful tuning of the improvement condition at each /3 is required 
so that physically the same condition is imposed. Chiral extrapolations and finite volume 
effects can also cause problems. 

Our values of c^, even on ensembles with high statistics, are rather imprecise after these 
considerations but improved errors would be possible with better smearings. Nevertheless 
it is clear that the smaller values of ca that we obtain, compared to those of the ALPHA 
collaboration, give improved scaling of the pion decay constant and the renormalised quark 
mass. 
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APPENDIX A: SIMULATION DETAILS 

The simulation details of the UKQCD quenched configurations are compiled in table |. 
Light quark propagators were generated using the Oia) improved Wilson action with the 
nonperturbative value of the improvement coefficient csw determined by the ALPHA col- 
laboration at /3 = 6.0 and 6.2 and the SCRI collaboration |0 at /3 = 5.93. In both cases 
the interpolating curves which were found to fit the nonperturbative determinations of cgw 
were used (rather than, for example, the numerical results at /3 = 6.0 and 6.2). The values 
of K were chosen in order to straddle the strange quark mass. The quark propagators were 
tied together to form mesons in both degenerate and non-degenerate k combinations. 

Gauge-invariant smearing was applied at the source and/or sink using either extended 
"fuzzed" spatial functions (denoted F) or Gaussian-like spatial functions using the 
Jacobi method (denoted S\ In the table, FL refers to a meson made up of a light quark 
propagator with a fuzzed source and local sink, combined with a LL antiquark propagator. 

A previous analysis of the data sets at /5 = 6.0 and 6.2 identified a small number 
of exceptional configurations: 3 on the small volume at /5 = 6.0, 2 for the large volume 
and 1 configuration at /5 = 5.93. These configurations were removed from the statistical 
ensemble. For more details see [§. The lattice spacing is set using tq = 0.5 fm and the 
interpolating curve determined by the ALPHA collaboration for ro/a [1^, from which we 



obtain ro/a = 4.741, 5.368 and 7.360 at f3 = 5.93, 6.0 and 6.2 respectively. Throughout the 
analysis we ignore all errors in tq. 
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APPENDIX B: EXTRACTING MASSES AND DECAY CONSTANTS 



We extract the mass and decay constant of the pseudoscalar meson from the two-point 
correlation functions, 



CpA, = (P^(0)A4(t)) (Bl) 

Cpp = (p\0)P{t)) (B2) 



where, 



oo -| 

CpA. = E Jt^T ((oI^V) (-1^410) e-^^^* - (0|pt|,\ ^n\A,\0) e'^^^^^-)) (B3) 



oo 



E^Ke-^^^/^sinh(M;5(t-|)) (B4) 



n=0 



Cpp = J2 Appe-*^?^/2 ^Qgh ( _ |) j (B5) 



n=0 



and Mp is the mass of the nth radial excitation of the pseudoscalar meson. At /3 = 6.0 and 
6.2, we performed a simultaneous correlated fit to 5 correlators, Cp^^, Cpp, Cpp, Cp^^ and 
Cpp, where FL etc refers to the smearing of the correlator. This combination was chosen 
in order to extract information on the excited states as well as the ground state masses 
and amplitudes. To increase statistics the correlators were averaged about the center of the 
lattice. 

Fits including the ground, first and second excited states were attempted. Obtaining 
these fits at /5 = 5.93 was straight forward. However, at /? = 6.0 the fits lie consistently 
below the central values of the 5 correlators, over any given fitting range, even with Q > 0.1. 
This is due to the fact that the small eigenvalues of the covariance matrix for the fits are not 
determined well enough with the statistics we have available. However, it is these eigenvalues 
which dominate the minimization of the x^- We chose to zero eigenvalues which where below 
a certain cut-off, CcutWmax, where Wmax is the largest eigenvalue and Ccut = 10"'^ for all k, 
combinations. This resulted in removing approximately 16 (22) eigenvalues for a 1 (2) state 
fit. 

The ground state masses obtained for the heaviest mesons at the 2 (3 values are shown 
in figure ^ as a function of the initial timeslice for the fit (tmin)- The results for Mpg when 
including only the ground state in the fit function were consistent with those including 
radial excitations, and for simplicity, we used these values in the final analysis shown in 



tables ||VI1| - |V111]| . The fitting ranges chosen were 11 — 16 at /5 = 5.93 and 13 — 24 at 
(3 = 6.0. The errors were generated using 1300 and 1000 bootstrap samples for jS = 5.93 and 
6.0 respectively. Our results are consistent, to within 2a, with the previous determinations 
of Mps in reference []8|,p!6[|, which were obtained using a subset of correlators used here. 

At (3 = 6.2, it was not possible to simultaneously fit to all 5 correlators and instead we 
used Cp^^, Cpp and Cpp to extract the ground state mass and decay constants and Cp^^, 
Cpp, Cp\^, Cpp to extract the ground state and first excited state contributions to the 
correlators being used to extract ca (see below). In both cases it was necessary to remove 
eigenvalues with Ccut = 10~^. For the 3 correlator case only a 1-state fit was successful, as 
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the smearing was optimised for the ground state and the contributions from excited states 



are small. The ground state mass as a function of t^m is shown in figure [TT| for the heaviest 
and lightest k combinations. For the latter, the mass falls off steadily with tmin) even though 
the Qs for the fits are reasonable. Comparing the fit with the correlators, we find only at 
tmin = 19 can we be confident that residual contributions from excited states are below the 
statistical errors. This is true for all n, except the heaviest. We choose 19 — 25 for all us. 
The final values for Mp are given in table 11^. Fitting to 4 correlators, including LL, 1, 2 
and 3 state fits could be performed. 

The masses and amplitudes obtained from the 2-state fits can be used to determine the 
contribution of the ground and first excited states to the correlators CpA4 and Cpp, which 
are used in the calculation of ca- The time dependence of this contribution is interesting 
to compare with the fitting range chosen to extract ca- We substitute the parameters from 
the 2-state fits with the fitting range 7—15, 6 — 24 and 8 — 24 at = 5.93, 6.0 and 6.2 
respectively, into equation and calculate the ratio of this for Umax = 1 and 2 with the 
correlators CpA^ (similarly for Cpp). The fractional contribution of the ground state, and 
the sum of the ground and first excited states, as a function of timeslice is shown in figures |^ 
for /5 = 5.93. Note that for a 2-state fit the excited state masses and amplitudes are likely 
to contain some contamination from higher excited states and results in the figures can only 
give a rough indication of the fractional contributions. 

We see that for the LL (FL) correlators the ground state dominates at around timeslice 
tgs = 12 (12), while the first excited state becomes the dominate radial excitation from 
approximately timeshce tist = 6 or 7 (5). Repeating this analysis at /3 = 6.0, we find tgs = 12 
for LL correlators and tist = 7. The smearing at this (3 was optimised for the ground state 
and hence it dominates much earlier, roughly timeslice 8 for the FL correlators; tist = 5 and 
the fractional contribution of the first excited state is very small compared to that for the 
LL. At (3 = 6.2, the smearing was optimised in a similar way and tgs = 10 for FL compared 
to 15 for LL correlators; tist ~ 8 and 6/7 for LL and FL correlators respectively. We see 
that the 1st excited state dominates at roughly the same physical time; tist = 6 at /5 = 5.93 
corresponds to t = 7 at /3 = 6.0 and ~ 9 at /5 = 6.2. We now concentrate on the ground 
state mass and amplitude only and drop the superscript in the following. 

The bare improved pseudoscalar decay constant, equation |l], 

fimp ^ j(0) ^ (gg) 

can be obtained from the amplitudes and masses: 



where 



or 



a/m = A^/<») (B8) 



A = sinh(aMp) (B9) 
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A = sinh(aMp) - -(2 sinh(aMp) - sinh(a2Mp)) (BIO) 
6 

depending on whether we have used the O(a^) or 0{a'^) definition of 84, respectively, to be 
consistent with the determination of ca- The values obtained for f^^^ and af^^^ extracted 
from the 1-state fits are given in tables ||VIl| - pC| . For (3 = 5.93 and 6.2 the values of 



are also given, where ca = —0.032(14) and —0.031(5) was used, respectively (obtained by 
chirally extrapolating ca - see next section). The statistical errors in ca are included by 
bootstrapping, into the error estimate of 

The bare improved PCAC mass is given by equation 0, 

^PCAC — ^PCAC + CaO'^pcAC- (-^11) 

For consistency throughout the analysis we extract f^^^cAc and fT^\^cAc using the masses 
and amplitudes extracted in the 1-state fits described above. 

m%AC = (B12) 

am?cAc = ^A^'^ (B13) 

where 

A(^) = 2[cosh(aMp) - 1] (B14) 



and 



A(2) = 2[cosh(aMp) - 1] - i[cosh(2aMp) - 4cosh(aMp) + 3] (B15) 

6 



for O(a^) or ©(a'^) temporal derivatives respectively. The results are detailed in tables [ VTl 



IX| . am^^^ is also given for P = 5.93 and 6.0. The values for m'^°^ and m^^^ are consistent 



to within 2a with those obtained by performing a constant fit to r(t) and s{t) directly. 



APPENDIX C: RENORMALISED DECAY CONSTANT AND QUARK MASS 

The renormalised decay constant is obtained from the combination 

r^ = {l + abAmg)ZAn''. (CI) 

The factor Za has been calculated nonperturbatively by the LANL P] and ALPHA 
groups (as well as by other groups using axial- Ward identities [|1^ and other methods [pJ|] ). 
Their results are consistent to within 2a at (3 = 6.0 and 6.2. The ALPHA collaboration 
investigated several (3 values in the range 6.0 to 24.0 and obtained the interpolating curve, 

_ l- 0.8496^g + 0.0610^0^ 
^ 1 - 0.7332^2 • I ) 
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In perturbation theory Za does not depend on ca to 1-loop. In addition, the perturbative 
value of Za is within a few percent of the nonperturbative value at /3 > 6.0 and hence we do 
not expect Za to depend significantly on ca when determined nonperturbatively over our 
range of /?s. Thus, we use the interpolating curve above for Za- For the errors in Za-, we use 
those from the direct simulations of the ALPHA group at /5 = 6.0 and 6.2 [Za = 0.7906(94) 
and 0.807(8) (2) respectively), and at /? = 5.93 we assign the same error as at /5 = 6.0. 
Although the interpolating curve is only valid in the range (3 > 6.0, we do not expect to 
incur a significant error by applying it at /5 = 5.93 since Za is not rapidly changing. 

The coefficient Ba has only been determined nonperturbatively at /5 = 6.0 and 6.2 by 
the LANL group They obtain Ba = 1.28(3) (4) and 1.32(3) (4) respectively compared to 
1.38(7) and 1.34(5) from 1-loop tadpole-improved perturbation theory where 

6a = — (1 + 0.8646a) (C3) 
Mo 

and we use a = ap(l/a) ^ and the fourth root of the plaquette for uq. We assign the 
error in the perturbative result to be lap(l/a). The perturbative and nonperturbative 
determinations of Ba are consistent and we use the perturbative results for our 3 (3 values. 
Finally, we take the quark mass appearing in equation |Cl] to be rnpcj^Q. The results for 
rof^'^"' at the reference mass r^Mpg = 3.0 are given in table Xlll , where we have applied the 



same values of Za and 6a for results obtained using different values of ca- 

Values for a renormalised quark mass are normally quoted jn_terms of the MS scheme at 
a particular reference scale, we choose 2 GeV. We calculate m^^'^ in two ways: nonperturba- 
tively, using the method suggested by the ALPHA collaboration, where one first calculates 
the renormalisation-group invariant quark mass 

M = ZMmp^Ac (C4) 



and then converts to the MS scheme at 2 GeV ||rT| using (up to) 4-loop perturbation 



theory. The ALPHA group have calculated Zm nonperturbatively I^Tf and obtained the 
interpolating curve 



1.752 + 0.321(/3 - 6) - 0.220(/5 - 6)^ (C5) 



for the range 6.0 < /3 < 6.5. The associated (/3 dependent) uncertainty in Zm is 1.1%. 
We ignore the additional, (3 independent error of 1.3% as we are only interested in scaling 
behaviour and not predictions in the continuum limit. The determination of Zm does depend 
on Ca through the extrapolation to zero quark mass (this limit is found using rnpcj^^ to 
determine Kc). The results of the next section will show that we obtain consistent results 
for Kc using our values of ca and those of the ALPHA collaboration and hence we do not 
expect a significant error from applying their values for Zm in our analysis. The values 
for rom^^"^(2 GeV) for r^Mpg = 3.0 obtained in this way are shown in table pClll|. The 



conversion factor from the renormalisation invariant mass to the MS scheme at 2 GeV is 
0.72076 at 4- loops [|rT|. We apply the same values of Zm for results obtained using different 

values of ca- 

We also calculated m*^'^(2 GeV) perturbatively using 
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m 



MS 



Za{1 + abAmg) 
Zp{fi){l + abpTUg) 



m 



imp 
PCAC 



where to 1-loop in tadpole-improved perturbation theory [ p2| , pU| ,P[ 

Za = no (1 - 0.416a) 
hp = —{l + 0.8763a) 



Zp(/i) 



Mo 
Mo 



1 + a(— ln(/ia)^ 
47r 



1.328) 



(C6) 

(C7) 
(C8) 

(C9) 



We assign errors to the perturbative resuhs of lap(l/a), as before. Using perturbative 
factors enables a meaningful comparison with the quark mass determined using (see next 
section), for which the renormalisation factor is only known perturbatively. The resulting 
values of m^^^ {2GeV) are given in table pClli. 



APPENDIX D: CHIRAL EXTRAPOLATIONS. 

The improvement scheme we have chosen defines ca in the limit of zero quark mass. We 
extrapolated our results for ca at finite quark mass using the fit function: 

CA{mqi,mq2) = Co + ci(mgi + mg2) + C2(mgi + mq2f + C3mgimg2, (Dl) 

where rriqi = Mp^i^ni), the pseudoscalar meson mass for degenerate quarks with 
We performed correlated fits for all tex, tgs combinations which gave reasonable Q values at 
finite quark mass. Higher order terms in the fit function were included successively, starting 
with a constant fit. Each fit function was applied to the set of CA{mqi,mq2) starting with 
the lightest m^i and mq2- 

The results are detailed in tables 0, [V| and [V^ for /? = 5.93, 6.0 and 6.2, respectively, 
where for each tex and tgs combination the fits which cover the largest quark mass range are 
given. In the case of a constant fit, ca{kc) is an average of the values at finite quark mass; 
we prefer to quote the value for the heaviest k value, which is consistent but has a larger 
error. 

At (3 = 5.93, we encountered difficulties performing the chiral extrapolations. Often the 
fit would lie above or below all the data points. In most cases, the problem was solved by 
fitting to only the degenerate k combinations, and hence reducing the correlations between 
data points. If this was not sufficient we also dropped the smallest eigenvalues from the 
covariance matrix, as described in the previous section. 

In some cases after eigenvalues were dropped there were not enough degrees of freedom 
remaining for the fit to be performed. Instead we performed a linear and quadratic (cs = 0) 
uncorrelated fit to the full quark mass range. These results are also given in the tables. For 
the comparison of ca{hc) as a function of t^x (see figure we take the result from the linear 
uncorrelated fit, unless there is a statistically significant difference between the two fits, in 
which case we use the quadratic fit. We also give in the tables the results of uncorrelated 
fits for comparison in the cases where we could not fit over the full quark mass range using a 
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correlated fit. The final value for ca at each /3 is chosen from the set of successful correlated 
fits only. 

Compared to the statistical errors in ca at finite fitting to only a limited number of 
data points gave rise to large errors in the chiral limit. At /5 = 6.0 and 6.2, we followed the 
same procedure. The correlations do not seem to be so severe, however, and it was possible 
to fit to the non-degenerate and degenerate k combinations and no eigenvalues needed to 
be set to zero. 

Using c^(kc); we obtain m^^^* perform a consistency check by chirally extrapolating 
rgMpg as a function of ^^o'^poAc! froin the definition of equation |^, and equations 
and |B13|, Mps and vupcac should vanish at the same point. In the first set of fits we 



used the same functional form as in equation ^T] with m^i = '"^pcac'I'^i)' ^^"^ used the 
same procedure as above. We repeated the analysis using ca as determined by the ALPHA 
collaboration. The results are given in table PCII| for /5 = 5.93 and ^ for /? = 6.0 and 6.2. 



At /3 = 5.93 we see that whether using a linear or quadratic fit function (with successively 
wider ranges in quark mass), there is a non-zero value for cq. Fits with cq forced to zero 
were unsuccessful. This effect was also seen at /5 = 6.0 and 6.2. We therefore tried to 
resolve the presence of quenched chiral logarithms. Following the analysis of the CP-PACS 
collaboration P3 we computed the quantities 



x = 2 + ^^^^i±^ln(^) ' (D3) 

- mg2 mgi 

which are related hy y = 1 + 6 ■ x + O(m^). Any significant deviation of y from 1 indicates a 
non-zero value for 6, the chiral log term. Figure |1^ shows that we observe a clear non-zero 
slope with 6 very roughly in the range 0.1 — 0.2, with y deviating from 1 by more than 3crs 
at /9 = 5.93. Previous estimates of 6 vary in the range of 0.08—0.12 from CP-PACS ||23|| , 



0.06 from Bernard et. al. |2^ and 0.065(13) from Bardeen et. al. [^. This analysis 



motivates us to include a log term in the fitting function in order to set Cq = 0.: 

mg2 - mgi mg2 

+C2{mgi + mg2f + Csmgimg2 (D4) 



The results of these fits are also shown in tables ^ and We present the best fits in 



figures |T2| using both our values for ca and those of the ALPHA collaboration. In figure (a) 
a quadratic fit (function 3 to the 8 lightest data points) to the /? = 5.93 data set is included to 
illustrate the non-zero intercept found when there are no chiral logs terms in the fit. Since 
rnpcj^Q is a bare mass there is no significance in the fact that there is better agreement 
between the results at different (3 values when our values for ca are used compared to those 
of the ALPHA collaboration. 

We also extrapolated rripcAc ^^^1 Mpg as function of 1/k in order to extract k^. The 
latter can be used to obtain independent of ca- There are no chiral logarithms expected 
for the PCAC mass and we use equation pi| for the extrapolations with cq = and = ihg 
defined as 



17 



m'(l + bmm'' 



q 

1 /I 



K 



1 



(D5) 
(D6) 



where Kc is an additional parameter in the fit. bm is known to 1-loop in perturbation 
theory Including tadpole-improvement 



-0.5 - 0.686a. 



(D7) 



We obtain bm = —0.730, —0.761 and —0.776 at /3 = 6.2, 6.0 and 5.93, respectively using 
ap{l/a). The results of the fits to TnpcAc shown in table PCI1| . To extract Kc from Mp 
we use equation with equation D3. The results of the fits are also given in table Xll . 
For comparison we present the results of fits without the chiral logs. 

We see a significant drop in the value of at /? = 5.93 and 6.0 from Mp^ when log 
terms are included. For 13 = 6.2 the results are consistent to within 2a. The values of Kc 
extracted from the PCAC mass and Mp should be consistent. We see that this is the case 
when comparing the PCAC mass results with those of the chiral log fits to the pseudoscalar 
meson mass. Disagreement between the two determinations of Hc when log terms are omitted 
for Mpg has been noted and discussed previously in reference and also in reference . 

Kc has been extracted previously by UKQCD from the simulations at /3 = 6.0 and 
6.2 in reference |0] and (5 = 5.93 in reference |T^. In those works uncorrelated linear 
extrapolations (without chiral logs) were performed to Mpg using all k combinations and 
with bjn calculated using boosted perturbation theory (the dependence of on the value 
of bm used was investigated and found to be small), = 0.135202(11), 0.135252t^^ and 
0.135815111 was obtained at /? = 5.93, 6.0 and 6.2, respectively. These values are consistent 
to within 2a with our results without chiral logs. 

The ALPHA collaboration has determined Kc from rripcAc at /5 = 6.0 and 6.2 [§: they 
find Kc = 0.135196(14) and 0.135795(13) respectively. These results are consistent, within 
2a, with those in table PCII| obtained using the ALPHA value for ca- 



APPENDIX E: RENORMALISED QUARK MASS FROM kc 

From the extrapolations of r^Mpg versus rorhq (equation P5|) we can extract the value 
of rhg which corresponds to r^Mpg = 3.0. This can be converted to a value for m^^ using 



7MS 



{afi)m. 



(El) 



The renormalisation factor is only known 1-loop in perturbation theory . Using tadpole- 
improvement. 



Zm{l/a) = — [1 + 1.002a] . 



(E2) 



The results for m {2GeV) are given in table XIII . 
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TABLES 



TABLE L Simulation details. 



p 


Volume 


^confs 


Csw 


smearing 


1^1 


L (fm) 


5.93 


16^ X 32 


684 


1.82 


LL, FL, FF 


0.1327, 0.1332, 0.1334, 0.1337,0.1339 


1.7 


6.0 


16^ X 48 


496 


1.77 


LL, FL, FF 


0.13344, 0.13417, 0.13455 


1.5 




32^ X 64 


70 




S L ^ S S 




3.0 


6.2 


243 X 48 


214 


1.61 


LL, FL, FF 


0.1346, 0.1371, 0.13745 


1.6 
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TABLE II. Details of the calculation of Cyi at /3 = 5.93 for the heaviest k value, k = 0.1327 
and FL and LL correlators. Standard, 0{a?) derivatives are used for the temporal derivatives. 



d = 0{a^) 





igs 




ritex)-ritgs) 


CA 


CA LANL 


Q 










Ar 
As 


fit ^ex to tgs 




FL 


Q 






n n770fi 1 'I 


-0 1 "^QO^ 




n 


A 




n 348f^'l 


n n354fQ'i 


-0 ^()9(9) 


-fl 101 c^^ 


n 






f) 181 (d) 


ni ^qfR^i 


-0 088('3') 


-0 08Ql'3') 


09 


6 




0.092(4) 


0.0069(6) 


-0.074(7) 


-0.078(6) 


0.04 


7 




0.052(3) 


0.0027(6) 


-0.051(11) 


-0.058(12) 


0.08 


8 




0.023(3) 


0.0012(5) 


-0.052(24) 


-0.065(33) 


0.09 


9 




0.014(3) 


0.0014(5) 


-0.104(48) 


-0.134(55) 


0.30 


3 


14 


0.552(6) 


0.078(1) 


-0.141(2) 


-0.130(2) 


0. 


4 




0.348(5) 


0.0365(9) 


-0.105(2) 


-0.101(3) 


0. 


5 




0.181(4) 


0.0170(7) 


-0.094(4) 


-0.090(3) 


0.04 


6 




0.092(3) 


0.0079(6) 


-0.086(7) 


-0.081(6) 


0.05 


7 




0.052(3) 


0.0037(6) 


-0.073(11) 


-0.068(11) 


0.05 


8 




0.022(3) 


0.0023(6) 


-0.101(28) 


-0.091(33) 


0.05 


9 




0.013(3) 


0.0025(6) 


-0.188(67) 


-0.162(57) 


0.39 


LL 


3 


12 


1.78(1) 


0.437(3) 


-0.246(1) 


-0.267(3) 


0. 


4 




0.758(8) 


0.134(2) 


-0.177(1) 


-0.177(3) 


0. 


5 




0.314(5) 


0.0416(8) 


-0.132(2) 


-0.133(3) 


0. 


6 




0.142(4) 


0.0141(7) 


-0.100(5) 


-0.105(5) 


0. 


7 




0.073(4) 


0.0048(7) 


-0.066(9) 


-0.073(8) 


0.09 


8 




0.032(3) 


0.0018(6) 


-0.058(17) 


-0.067(20) 


0.05 


9 




0.018(3) 


0.0018(G) 


-0.102(30) 


-0.125(39) 


0.32 


3 


14 


1.78(1) 


0.438(3) 


-0.246(1) 


-0.266(3) 


0. 


4 




0.757(8) 


0.135(2) 


-0.178(1) 


-0.176(3) 


0. 


5 




0.314(5) 


0.0427(9) 


-0.136(3) 


-0.132(3) 


0. 


6 




0.141(4) 


0.0152(7) 


-0.108(5) 


-0.106(4) 


0.0 


7 




0.072(4) 


0.0059(6) 


-0.082(9) 


-0.079(8) 


0.07 


8 




0.031(4) 


0.0030(6) 


-0.096(20) 


-0.082(19) 


0.04 


9 




0.017(3) 


0.0029(6) 


-0.174(49) 


-0.145(40) 


0.33 
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TABLE III. The same as for table || using improved (0(a^)) lattice derivatives. 





tgs 




ritex)-r{tgs) 


CA 


CA LANL 


Q 










Ar 
As 


fit t^x to tgg 




FL 


4 


12 


0.296(5) 


0.0103(10) 


-0.035(3) 


-0.036(3) 


0.05 


5 




0.149(4) 


0.0048(7) 


-0.032(5) 


-0.033(4) 


0.04 


6 




0.077(4) 


0.0017(7) 


-0.022(9) 


-0.024(9) 


0.04 


7 




0.045(4) 


-0.00002(71) 


-h0.0004(16) 


-0.012(16) 


0.03 


8 




0.018(3) 


-0.00009(61) 


-F0.005(35) 


-0.035(817) 


0.02 


4 


14 


0.295(5) 


0.0116(10) 


-0.039(3) 


-0.037(2) 


0.07 


5 




0.148(4) 


0.0061(8) 


-0.041(5) 


-0.035(4) 


0.05 


6 




0.076(3) 


0.0030(7) 


-0.040(10) 


-0.029(9) 


0.04 


7 




0.045(4) 


0.0013(6) 


-0.030(15) 


-0.025(15) 


0.02 


8 




0.017(3) 


0.0012(6) 


-0.072(45) 


-0.066(264) 


0.03 


LL 


4 


12 


0.483(6) 


-0.0276(12) 


-h0.057(3) 


-h0.034(2) 


0. 


5 




0.227(5) 


0.0010(8) 


-0.004(4) 


-0.007(3) 


0.01 


6 




0.112(4) 


0.0020(8) 


-0.018(7) 


-0.020(6) 


0.03 


7 




0.062(4) 


0.0001(8) 


-0.002(13) 


-0.012(12) 


0.02 


8 




0.025(4) 


-0.0001(6) 


-h0.003(26) 


-0.026(33) 


0.02 


4 


14 


0.481(6) 


-0.0262(12) 


-h0.054(3) 


-h0.032(2) 


0. 


5 




0.225(5) 


0.0024(9) 


-0.011(4) 


-0.010(3) 


0.0 


6 




0.110(4) 


0.0034(8) 


-0.031(7) 


-0.024(6) 


0.03 


7 




0.060(4) 


0.0015(7) 


-0.025(12) 


-0.022(11) 


0.02 


8 




0.024(4) 


0.0013(7) 


0.056(31) 


-0.046(30) 


0.01 
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TABLE IV. ca in the chiral limit for /3 = 5.93 using both O(a^) and O(a^) derivatives. Columns 
1—7 give details of correlated chiral extrapolations of ca determined from mesons with degenerate 
quark masses. Included are the function used (in terms of the index of the highest order coefficient 
included - see equation the number of eigenvalues dropped from the covariance matrix, Tid^op 
and the number of data points used in the fit, n^. In some cases a linear and quadratic uncorrelated 
fit (function 1. and 2., respectively) was performed. The results for ca(kc) for the uncorrelated 
fits are given in columns 8—9. 



tex 


tgs 


func. 






Q 


^drop 


uncor func 1. 


uncor func 2. 


FL 


5 


12 












-0.068(6) 


-0.079(10) 


6 




1. 


4 


-0.048(10) 


0.72 





-0.046(15) 


-0.037(27) 


7 




2. 


4 


-0.002(29) 


0.78 





+0.006(30) 


+0.036(56) 


5 


14 












-0.069(6) 


-0.079(10) 


6 




2. 


4 


-0.045(21) 


0.68 





-0.049(15) 


-0.040(28) 


7 




2. 


4 


+0.023(49) 


0.73 





-0.006(31) 


+0.033(61) 


LL 


7 


12 


2. 


4 


-0.006(30) 


0.74 





-0.019(22) 


+0.002(40) 


8 




1. 


4 


+0.011(50) 


0.98 





+0.016(61) 


+0.046(149) 


7 


14 


2. 


4 


-0.030(23) 


0.45 





-0.029(21) 


+0.000(39) 


8 




1. 


4 


+0.001(64) 


0.98 





-0.002(68) 


+0.068(295) 




FL 


4 


12 


1. 


5 


-0.050(3) 


0.05 









5 














-0.051(8) 


-0.063(13) 


6 














-0.018(21) 


+0.000(38) 


4 


14 


1. 


5 


-0.050(3) 


0.2 


1 






5 














-0.052(8) 


-0.064(13) 


6 




2. 


5 


-0.022(19) 


0.78 









LL 


5 


12 


1. 


4 


-0.039(4) 


0.01 





-0.041(6) 


-0.047(10) 


6 




2. 


5 


-0.028(13) 


0.85 









6 


14 


2. 


5 


-0.032(14) 


0.70 
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TABLE V. CA in the chiral limit for /? = 6.0. The details are the same as in table No 
eigenvalues were dropped from the covariance matrix. In the case of a constant fit the (*) indicate 
the value of ca from the heaviest n value is taken rather than the value from the fit. 





tgs 


form 






Q 


u lin 


u quad 


FL 


4 


12 


1. 


4 


-0.124(58) 


0.94 


-0.095(52) 


-0.154(120) 


4 


16 


1. 


4 


-0.145(59) 


0.95 


-0.113(54) 


-0.188(153) 


LL 


6 


12 


2. 


4 


-0.084(53) 


0.69 


-0.056(17) 


-0.069(33) 


7 




1. 


5 


-0.106(48) 


0.97 


-0.114(52) 


-0.166(172) 


6 


16 


1. 


4 


-0.065(12) 


0.99 


-0.063(16) 


-0.067(30) 


7 




3. 


6 


-0.147(66) 


0.91 






SL 


4 


12 


1. 


6 


+0.002(12) 


0.1 






5 







6 


+0.008(16) 
+0.006(18)* 


0.17 






4 


16 


3. 


6 


-0.008(19) 


0.92 






5 




1. 


5 


+0.017(25) 


0.74 


+0.014(40) 


+0.101(77) 



TABLE VL CA in the chiral limit for /? = 6.2. The details are the same as in table IV. No 
eigenvalues were dropped from the covariance matrix. In the case of a constant fit the (*) indicate 
the value of ca from the heaviest k value is taken rather than the value from the fit. 







tgs 


form 






Ca(Kc) 


q 


FL 


4 


16 


1. 




6 


-0.021(11) 


0.65 


5 




1. 




6 


-0.048(35) 


1.0 


LL 


6 


16 


3. 




6 


-0.043(7) 


0.75 


7 




0. 




6 


-0.029(4) 
-0.031(5)* 


0.76 


8 









6 


-0.028(7) 
-0.029(9)* 


0.89 
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TABLE VII. The pseudoscalar meson masses, decay constants and ampcAC extracted at 
(3 = 5.93. a?"rn'plj and a? f^^^ were extracted using O(a^) derivatives. The error in ca is in- 
cluded in the results for a'rnpcAC ^"^^ a/*™^. 







^2 


aMps 


(0) 


2 (1) 

a rnp(jj^(j 


imp 




a/(0) 




^pmp 


0.1327 





1327 


0.4948(10) 


0.0727(1) 


0.1223(5) 


0.0688(18) 





131(1) 


0.219(2) 


0.124(3 


0.1332 





1327 


0.4684(10) 


0.0655(1) 


0.1096(5) 


0.0619(16) 





128(1) 


0.214(2) 


0.121(3 


0.1334 





1327 


0.4575(10) 


0.0626(1) 


0.1046(5) 


0.0592(15) 





127(1) 


0.212(2) 


0.120(3 


0.1332 





1332 


0.4409(11) 


0.0583(1) 


0.0972(5) 


0.0551(14) 





126(1) 


0.209(2) 


0.119(3 


0.1334 





1332 


0.4296(11) 


0.0553(1) 


0.0922(5) 


0.0524(13) 





125(1) 


0.207(2) 


0.118(3 


0.1334 





1334 


0.4179(12) 


0.0524(1) 


0.0873(5) 


0.0496(13) 





123(1) 


0.205(2) 


0.117(3 


0.1337 





1337 


0.3814(13) 


0.0436(1) 


0.0727(5) 


0.0413(11) 





120(1) 


0.200(2) 


0.114(3 


0.1339 





1337 


0.3686(13) 


0.0407(1) 


0.0679(5) 


0.0385(10) 





119(1) 


0.198(2) 


0.112(3 


0.1339 





1339 


0.3553(14) 


0.0377(1) 


0.0631(5) 


0.0357(9) 





118(2) 


0.197(2) 


0.111(3 



TABLE VIII. Same as in table ^ for /3 = 6.0 







aMps 


(0) 


2 (1) 

a mp(j^(j 






0.13344 


0.13344 


0.3979(10) 


0.0532(4) 


0.0791(4) 


0.111(1) 


0.164(2) 


0.13417 


0.13344 


0.3555(12) 


0.0425(4) 


0.0632(4) 


0.107(1) 


0.159(2) 


0.13455 


0.13344 


0.3317(14) 


0.0366(3) 


0.0550(4) 


0.105(1) 


0.157(2) 


0.13417 


0.13417 


0.3078(13) 


0.0317(3) 


0.0474(4) 


0.103(1) 


0.153(2) 


0.13455 


0.13417 


0.2801(15) 


0.0258(2) 


0.0392(4) 


0.100(1) 


0.152(2) 


0.13455 


0.13455 


0.2493(18) 


0.0201(2) 


0.0311(4) 


0.098(1) 


0.151(3) 



TABLE IX. The same as in table for (3 = 6.2 





H2 


aMps 


(0) 


2 (1) 

a rup^^^ 


imp 
^^PCAC 


a/(0) 




^pmp 


0.1346 


0.1346 


0.2798(17) 


0.0363(1) 


0.0391(5) 


0.0351(2) 


0.079(1) 


0.085(2) 


0.076(1) 


0.1351 


0.1346 


0.2484(19) 


0.0289(1) 


0.0309(5) 


0.0279(2) 


0.076(1) 


0.081(2) 


0.073(1) 


0.1353 


0.1346 


0.2351(21) 


0.0259(1) 


0.0276(5) 


0.0250(2) 


0.074(1) 


0.079(2) 


0.072(1) 


0.1351 


0.1351 


0.2144(22) 


0.0215(1) 


0.0230(5) 


0.0208(2) 


0.073(1) 


0.078(2) 


0.070(1) 


0.1353 


0.1351 


0.1997(23) 


0.0185(1) 


0.0199(5) 


0.0178(1) 


0.071(1) 


0.077(2) 


0.069(1) 


0.1353 


0.1353 


0.1834(26) 


0.0155(1) 


0.0168(5) 


0.0150(1) 


0.070(1) 


0.076(3) 


0.068(1) 
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TABLE X. Parameters from r'^M'pg vs rQin^j^^ at P = 5.93. The function L2 etc refers to 



the fit function in equation D4, 71^ to the number of data points in the fit and Udrop to the number 
of eigenvalues omitted from the covariance matrix. Where the errors on the coefficients are not 
symmetric, assymetric errors are quoted. 



fit 




Q 


'^drop 


Co 


Cl 


Clog 


C2 


C3 










CA 


= NP 








1 


3 


0.19 





0.08(3) 


8.2(2) 








2 


5 


0.24 





0.23(1) 


7.3(2) 




1.2(1) 




3 


8 


0.83 


2 


0.27(3) 


7.1(2) 




1.6(2) 


-1-1(1) 


L2 


5 


0.14 







7.5(2) 


1.0(2) 


2.3(3) 






9 


0.20 


2 




7.5(1) 


1.0(1) 


2.4(i) 












CA = 


ALPHA 








1 


3 


0.19 





0.17(3) 


8.9(6) 








2 


5 


0.27 





0.24(4) 


8.4(1) 




0.8(1) 




3 


8 


0.59 


2 


0.27(1) 


8.3(1) 




0.69(13) 


1-2(1) 


L2 


4 


0.20 







8.2(2) 


1.6(2) 


3.2(t) 




L3 


7 


0.27 


1 




8.4(1) 


1.3(2) 


1.8(2) 


3.4(5) 



TABLE XL Parameters from r^Mpg vs tqitiJ^qj^q at /? = 6.0 and 6.2. This time only for the 
best fits. The details are the same as in table 



fit 




Q 


^drop 


Co 


Cl 


Clog 


C2 


C3 










/3 = 6.0, CA = 


ALPHA 








3 


6 


0.76 





0.29(7) 


7.7(4) 




1.9(1.1) 


-3.2(1:3) 


L3 


6 


0.39 







7.7(4) 


1.7(4) 


3.7(}:i) 


-1.5(2.0) 










13 = 6.2, CA 


= NP 








3 


6 


0.14 





0.27(i) 


6.6(1) 




2.1(4) 




L3 


6 


0.58 







6.5(1) 


1.7(t) 


4.o(?:?) 


i.4(i:§) 










(3 = 6.2, CA = 


ALPHA 








2 


6 


0.14 





0.28(i) 


6.6(1) 




2-2(1) 




L2 


6 


0.67 







6.6(1) 


1.7(t) 


4.4(?:8) 
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TABLE XII. He determined from Mps and 


''^PCAC- 


The parameters 


"fit" and are the same 


as 


in table |Xl 
















CA = NP 


imp 
'"^PCAC 

CA = NP{ALPHA) 








= 5.93 






fit 


2 


L2 




2 


1 




4 


7 




4 


4 


Q 


0.39 


0.42 




0.41 


0.60 




0.135252(11) 


0.135089(1^) 




0.135126(11) 


0.135099(1) 






P 


= 6.0 






fit 


3 


L2 




- 


1 




6 


6 




- 


6 


Q 


0.77 


0.28 




- 


0.33 


Kc 


0.135291(if) 


0.135175(1^) 




- 


0.135185(^°) 






P 


= 6.2 






fit 


1 


L2 




1 


1 




5 


6 




6 


6 


Q 


0.35 


0.36 




0.34 


0.31 




0.135820(16) 


0.135786(i|) 




0.135816(4) 


0.135816(4) 



TABLE XIII. Decay constants and renormalised quark mass at rlMls = 3.0: (a) usmg CA 
from this paper, (b) using ca as determined by the ALPHA collaboration. The two errors given 
are statistical and systematic. In the case of (a) the statistical errors include those from ca, while 
for (b) the statistical errors from ca are much smaller and are not included. The systematic er- 
rors are as follows. For rom'^Qj^Q{2GeV) from Zm the errors correspond to a 1.1% error in Zm- 
For rQmJp(jj^fj{2GeV) from Za/Zp we assume a lap(l/a) error in Za/Zp (the mass corrections 
{bA — bp)amq are ignored as they are negligible). Similarly for Z^ used for rQm'^f {2GeV) . For 
^^yren ^j^g systematic error is dominated by the error in Za obtained by the ALPHA collabora- 



tion [Tj. 





6.2 


6.0 


5.93 




romrcAcC^G^V) (a) 


0.240(2)(3) 




0.227(6)(2) 


(b) 


0.238(1)(3) 


0.212(2)(2) 




Za/Zp 


rom^bAcC^GeV) (a) 


0.234(2)(11) 




0.219(6)(15) 


(b) 


0.232(1)(11) 


0.207(2)(13) 




rom'^f{2GeV) 


0.239(7)(13) 


0.232(4)(16) 


0.229(3)(17) 


ror^ (a) 


0.443(7) (4) 




0.440(12)(5) 


(b) 


0.439(6) (4) 


0.402(4)(5) 
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♦ ALPHA 

♦ LANL-2 

♦ LANLJ 
■ this paper 

^M|f TAD-PERT 



FIG. 1. Ca extracted by various groups. The 1-loop tadpole- improved values for ca were cal- 
culated using ap{\/a). Also shown as a black line is the Fade expansion of the ALFHA results 
The ALPHA and LANL results at /5 = 6.2 are offset for clarity. LANL-3 and LANL-2 refer to 
the results of the LANL group using standard O(a^) derivatives (equation 15) and modified O(a^) 
derivatives respectively. 
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FIG. 2. m{t) as a function of timeslice for (3 = 5.93, = 0.1327 and LL smearing. The dotted 
line indicates a fit to a constant. 
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FIG. 3. The fractional contribution of the ground state (squares) and the sum of the ground 
state and first excited state (diamonds) to the Cpa4 and Cpp, LL and FL correlators for (3 = 5.93 
and AC = 0.1327. 
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0.05 



0.00 



-0.05 



-0.10 



□ LL 5 = 0(a*) 
K LL 6 = 0(a^) 
O FL 5 = 0(a^)- 
o FL 6 = 0(a^) 
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FIG. 4. The values of ca obtained using the LANL method at /3 = 5.93 and ki = 0.1327 as a 
function of tex used in the fit. S indicates the type of lattice derivatives used. 
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FIG. 5. m{t) = r{t), s{t) and niimpit) = r{t) + acAs(t) as a function of timeslice at /3 = 5.93 and 
Ki = 0.1327 for LL correlators. The dotted lines indicate the values for 2ampcAC obtained from 
fitting the results to a constant The dashed lines indicate the value of 2ani^^^ obtained from the 
fit for CA, where the fitting range 7—12 was used to extract ca for O(a^) derivatives (bursts) and 
6 — 12 for O(a^) (squares). These data points also include the statistical error of ca- 
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0.05 



0.00 



n LL 6=0{a'^jss 

w LL (5 = 0(a'^) 

o FL (5 = 0(a*) 

o FL (5=0(a^) 
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FIG. 6. The values of ca obtained from chiral extrapolation in aMpg sX (5 = 5.93 as a function 
of tea; used in the LANL fit. The dashed points indicate that an uncorrelated chiral extrapolation 
was performed. 
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FIG. 7. CA as a function of aMpg at /? = 5.93 for LL correlators. The fitting range 6 — 12 was 
used for O(a^) derivatives (squares) and 7—12 for O(a^) derivatives (diamonds). 
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,= 16 



0.0 
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10 



(c) (d) 
FIG. 8. (a) CA as a function of tex for (3 = 6.0 and ki = 0.13344 with tgg = 12. O(a^) improved 
derivatives have been used. We found that the Qs for the LANL fits were higher than those 
obtained at /? = 5.93. (b) the same as in (a) but with tgg = 16. (c) ca in the chiral Umit from 



I'QS 



12. (d) ca{kc) from tgg = 16. 
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FIG. 9. (a) CA as a function of tex for /? = 6.2 and k = 0.13460 with tgg = 16. O(a^) improved 
derivatives have been used, (b) ca in the chiral hmit. 
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FIG. 10. The ground state masses of the pseudoscalar meson as a function of tmin for the 
heaviest mesons at /3 = 5.93 and 6.0. Cp^^, Cpp, Cpp, Cp\^ and Cpp were fitted simultaneously, 
using a fit function including ground (n = 1) and radially excited states (n = 2,3). The solid 
data points indicate fits for which Q > 0.1. For the dashed results Q > 0.01 but < 0.1. The 
star indicates the final value of tmin chosen. We found this fitting range to be adequate for all k 
combinations. 
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FIG. 11. The ground state masses of the pseudoscalar meson as a function of tmin for the 
heaviest and lightest meson at /3 = 6.2. Cp^^, Cpp and Cpp were fitted simultaneously, using a 
fit function including the ground state. The solid data points indicate fits for which Q > 0.1. For 
the dashed results Q > 0.01 but < 0.1. The star indicates the final value of tmin chosen. We used 
this fitting range for all k combinations. 
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i rorupcj^^Q. (a) using ca as determined in this paper and (b) as determined 
by the ALPHA collaboration, (c) shows the quantity y as a function of x as defined in equations 
and D3 for all 3 (3 values, (d) presents the results for /3 = 5.93 only. Note that in (c) and (d) 



FIG. 12. rgM|, 



''^PCAC been used. 
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FIG. 13. fQMpg and tqui^qj^^q vs rorriK- is calculated using Kc obtained from the fits 
shown for each (3 value. In (b) our results for ca are used and in (c) the ca determined by the 
ALPHA collaboration. 
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FIG. 14. The renormalised pseudoscalar decay constant, (a) r^f^^ vs rQMpg using the 
values of ca calculated in this paper. The errors are statistical, and include the statistical error 
of Ca- (b) as in (a) using ca from the ALPHA collaboration. Here the errors on ca have been 
neglected, (c) The scaling of tq/''^" for the reference mass {r^MpsY = 3.0. The dashed line 
indicates a simultaneous fit to the data sets obtained using our ca values and those from the 
ALPHA collaboration. The errors include the uncertainties in the coefficients bA and Za- The 
results at /3 = 6.2 have been offset slightly for clarity. 
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FIG. 15. Scaling of rQm^^{2GeV) for rQMpg = 3.0 (a) obtaining the renormalised quark 
mass from the bare PCAC quark mass via the renormalisation-group invariant mass (see equa- 
tion |C^ ) and using our values of ca and those of the ALPHA collaboration. The errors include the 
uncertainty in ca (for our ca values only) and Zm ■ The dashed line indicates a simultaneous fit to 
both data sets, (b) As in (a) but calculating m'^^ directly using the 1-loop perturbative results. 
This is compared and simultaneously fitted to the 1-loop results for the quark mass extracted using 
Kc, which is independent of ca- The errors include estimates of the perturbative uncertainty in 
Za/Zp. In both plots the data points have been offset for clarity. 
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